Non-linear self-consistent response of graphene in time domain 
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We study the non-linear electromagnetic response of graphene taking into account the self- 
consistent-field effects. Response of the system to a strong pulse excitation is calculated. It is 
shown that radiative decay in graphene differs from that of conventional two-dimensional electron 
systems both quantitatively and qualitatively. Possible applications of the predicted effects for 
generation of terahertz radiation are discussed. 

PACS numbers: 78.67.-n, 73.50.Fq, 81.05.Uw 

Recently discovered [1101 '^sw two-dimensional (2D) material - graphene - attracted much attention in the past three 
years |3| • Relativistic-Hke energy spectrum and the vanishing effective mass of charge carriers in graphene result 
in unconventional and interesting transport and electrod yna mic properties, as has been demonstrated in numerous 
experimental and theoretical studies 0,11, 0, B S E [IE 111 E, 111 [ll, H E H H 

3[2i,Sl3H. 

As a direct consequence of the "relativistic" spectrum, electromagnetic response of graphene was predicted [sij to 
be strongly nonlinear: Irradiation of a graphene layer by electromagnetic wave with the frequency should lead to 
the higher harmonics generation at the frequencies mfJ, m = 3, 5, 7, . . ., with the higher harmonics amplitudes falling 
down very slowly, as 1/|to|. 

In general, considering response of charge carriers to the external electromagnetic radiation, one should take into 
account self-consistent-field effects [s^l- The external time-dependent electric field of the wave Yf^^'it) induces in the 
2D gas, lying in the plane z = 0, the electric current ]{t). According to the Maxwell equations, this time-dependent 
current produces, in its turn, a secondary (induced) electric field 



Erio(i) 



-2^j(t)/c, 



(1) 



which is added to the external field and acts back on the electrons. It is the total self-consistent electric field 
E*°*Q(t) = E|^*g(<) -I- W^J^Q{t) that should be written in the equations of motion for the current j(i), instead of the 
external field ^%%{t). 

In conventional 2D systems with the parabolic electron dispersion, taking into account the self-consistcnt-field 
effects leads to the following equation for the current (j = — en^v, dv/dt — — eE*°*Q/m*) 



dt 



(2) 



where 



27rrise^/m*c is the radiative decay rate [SJ], Ug is the 2D electron density, m* is the electron effective 



mass, V is the hydrodynamic velocity of electrons, and the subscript par reminds that Eq. ^ refers to the conventional 
2D electrons with the parabolic energy dispersion. In ^ we have ignored the scattering due to impurities and phonons 
[the corresponding term 7j(t) can be added to the left-hand-side of ^] . In high-electron mobility samples the radiative 
decay rate Tpar substantially exceeds the impurities/phonons scattering rate 7, Tpar ^ 7- In finite-size 2D electron- 
gas samples the plasmon-, the cyclotron (in a magnetic field), and/or the magnetoplasmon resonances can be observed; 
the linewidth of these resonances in high-electron- mobility samples is mainly determined by Tpar 1^, 34 1 . 

The electromagnetic response of graphene is strongly nonlinear [3l[. In Ref. [sH, however, it was calculated 
ignoring the self-consistent-field effects. The goal of this Letter is to derive the corresponding equations of motion for 
the current j and to calculate response of the system to a pulse excitation E'^^*(i) = EoT5{t), where Eq and r are the 
amplitude and the duration of the pulse. The self-consistent non-linear response of graphene to a harmonic external 
electric field will be discussed elsewhere 



Like in [3l|, we consider a 2D electron gas with the "relativistic" massless energy spectrum ep± = ±V 



under the action of an external time-dependent electric field E'^^*(i). Assuming that the Fermi energy > lies 
in the electron band and the temperature is small, T <^ ep, we describe the graphene response using the kinetic 
Boltzmann equation in the coUisionless approximation, 



dt 



Vp/p(i)eE*°io(t) = 0. 



(3) 



2 



Here the sign + (— ) corresponds to the electron (hole) band, V « 10® cm/s is the effective "velocity of light" in 
graphene, and fp+{t) = fp{t) is the momentuni distribution function of electrons (from now on we omit the sign + 
for brevity). Equation ^ has the exact solution 

/p(0 =-^o(p-Po(t)), (4) 
where ^o(p) is the Fermi-Dirac function, and Po{t) resolves the classical single particle equation of motion 

^ = -eKMt). (5) 
The electric current j(i) = —egsgvS^^ vfp(t) is written in terms of po(i) as 



jW = J dp^dpy^To (p - po(t)) 



(6) 



where gs = .9t, = 2 are the spin and valley degeneracies in graphene, p — yPx '^Pyi ^ the sample area. 
Combining Eqs. (O, ([6]) and ([T]) we get the following equation of motion for the momentum po(t): 

dY>o{t) e^gsgvV f p -^ext /7^ 

+ ^ ^2 / dp^dpy-To (p - po(i)) = -eEz=S)- (7) 



dt 27:1%' 



c 



Equations ([7]) and ^ describe the non-linear self-consistent response of graphene to an arbitrary external electric 
field E|5*Q(i) in the considered approximations. After the non- linear equation ([7]) is resolved with respect to the 
momentum po{t), the current j(t) can be found from ([6]). 

Now assume that the temperature is zero and that the external electric field E*^^*(t) is directed along the x-axis. 
Then po = {pq,0) and the dimensionless x-component of the momentum P(t) = Po(t)/pp is determined by the 
equation 



dP{t) 



e 



dt pf 
where pp = £f/V is the Fermi momentum, 

-2 r/' 



rG(p) = i^r(^ = o,i), (8) 



G{P) = Vl + P^- cos 0de {y^l + Qcos9 - y^l-Qcose) , Q=- ^ < I, 

n Jq \ / 1 + 



(9) 



and 



gsgve^2eF .^e^ 



The function G{P) determines the electric current: 

3 



G(P). (11) 



In the linear-response regime, when the external electric field is so small, that |P| = \po{t)\/pF ^ 1; the function 
G(P) is linear, G(P) « P at |P| < 1, and Eq. ® gives 

^ + rP(t) = -—Sriz = 0, t). (12) 
dt Pf 

Equation (fT2|) is similar to ([2]) , and one sees that the frequency F has the physical meaning of the radiative decay rate 
in graphene in the linear regime. In contrast to Tpar, F is proportional to the square-root of the charge carrier density. 
For experimentally relevant Us the value of F lies in the subterahertz range, F/27r(TIIz)« OAly/usilO^'^ /cm^). 

Now consider response of graphene to a pulse external electric field E^^*{z = 0,t) = EorS^t), not imposing any 
restriction on the field amplitude Eq. Eq. ^ is reduced to the homogeneous differential equation 



3 



with the boundary condition 



f]P 

_+rG(P) = 0, t>0, (13) 



P(i^+0)^P„^-^. (14) 

PF 



The problem (fT3| -(fT4 |) can be solved analytically. Expanding the integral in ^ in powers of Q = 2P/(1 + P^) < 1, 
we get for the function G{P): 



G{P) 



Vp^TT 



3 f 2P V f 2P 



32 V -P^ + 1 / 1024 V -P^ + 1 



(15) 



The expansion ([T5|) is valid both at P ^ 1 and P ^ 1; even at P w 1 Eq. (fT5|) is very accurate: the third term is 
square brackets gives a correction smaller than 3.5%. Integrating now Eq. (|13p with the boundary condition (|14p we 
finally get the following implicit relation between the dimensionless momentum P{t) and the time t: 



jj;—- , 1 + 3/8 13/96 13/160 

V T -L 111 p, ^ ^)l/2 (p/2 _^ 1)3/2 (p/2^ 1)5/2 




Figured^ shows the dependence P(i), given by Eq. (fT6|l . at different values of the electric field parameter Pq = 
eEgT/pp. If the external field is small, Pq <C 1, the system relaxes after the pulse excitation exponentially, similar to 
the conventional 2D systems with the parabolic dispersion, 

P(t) -Poexp(-rt), Po<l. (17) 

The characteristic decay time is determined in this case by the inverse radiation decay rate PH)) . If the external field 
is strong, Pq 3> 1, the momentum of the system decays linearly in time, 

p(i) = Po-rt, Po>i. (18) 

The linear dependence p8)) remains valid until P{t) reduces down to P{t) ~ 1 (until t ~ Pq/F); after that P{t) decays 
exponentially like in pT|) . 

The current j(t) in the strong excitation regime Pq 3> 1 is equal to its highest possible value jmax = {—e)nsV 
and time-independent a,t t ^ fo/r, and then quickly (exponentially) decays (at the time scale ~ F"-'^) down to zero. 
Figure [T]3. 

The fact, that after the pulse excitation electrons in graphene move with a constant velocity V ~ 10^ cm/s for quite 
a long time ~ Po/F, may have interesting applications. In a finite-size graphene sample such excited electrons will 
be reflected by the boundaries and oscillate in the sample with the typical frequency ^ V/L, lying in the terahertz 
range, if the sample dimensions L ^ 1 fim. As at Pq S> 1 the time Pq/F is much longer than the oscillation period, 
this may lead to a coherent terahertz radiation from graphene excited by a strong pulse electric field. 

To conclude, we have derived equations describing the non-linear self-consistent response of graphene electrons to 
an external time-dependent electric field and calculated response of graphene to a pulse excitation. We have shown 
that at low excitation strengths, the system responds exponentially, like conventional 2D electron layers, with however 
different characteristic decay rate. The radiative decay rate F in graphene (fTO|) is proportional to the square root 
of the electron density, in contrast to the normal 2D systems, where such a dependence is linear. At the strong 
excitation strengths, the average momentum of graphene electrons decays linearly, in contrast to the exponential 
decay in conventional systems, with the average current remaining constant during the time ~ Po/F. The predicted 
effects may be used for terahertz generation. 

This work has been partly supported by the Swedish Research Council and INTAS. 
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FIG. 1: (Color online) The time dependence of (a) the momentum P{t)/Po and (b) the electric current j{t)/{—e)nsV , at a 
pulse excitation of graphene. Different curves correspond to different pulse amplitudes Po ~ {—e)EoT /pp. 
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